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Abstract 

A natural procedure is introduced to replace the traditional, per- 
turbatively generated counter terms to yield a formulation of covari- 
ant, self- interacting, nonrenormalizable scalar quantum field theories 
that has the added virtue of exhibiting a divergence-free perturbation 
analysis. To achieve this desirable goal it is necessary to reexamine 
the meaning of the free theory about which such a perturbation takes 
place. 

1 Introduction 

Nonrenormalizable quantum field theories, such as ip^ models, with a space- 
time dimension n > 5, need nontrivial counter terms for otherwise they lead 
to (generalized) free theories as shown by Aizenman [1] and Frohlich [2]. A 
free quantum theory has a trivial classical limit and so it can not correspond 
to the correct quantization of the original nontrivial classical theory. Regu- 
larized perturbation theory suggests an unending series of distinct and ever 
more singular counter terms which cannot be considered an acceptable so- 
lution. This situation suggests that we look elsewhere for suitable counter 
terms, and this paper addresses one such search. 

In Sec. 2 we outline a Euclidean space lattice model for ip^ models for 
n > 5, which includes an unconventional counter term along with the tradi- 
tional terms expected in such a lattice formulation. In Sec. 3 we present the 
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background for choosing this form for the model and explain the rationale for 
choosing the atypical counter term. Section 4 shows that general correlation 
functions can be suitably bounded by correlation functions at a sharp time 
as determined by the ground-state distribution. In Sec. 5 we take up the 
question of the continuum limit and study mass renormalization, coupling 
constant renormalization, and field strength renormalization. Here we show, 
thanks to the properties of our chosen counter term, that a suitable pertur- 
bation theory of the quartic interaction is divergence free. Importantly, this 
perturbation theory is not about the usual free theory but about a pseud- 
ofree theory, which is a model that contains the atypical counter term but 
docs not include the quartic interaction. The reason for this divergence-free 
character is related to a simple idea already illustrated by idealized lattice 
space integrals presented in Sec. 3. Finally, Sec. 6 offers additional discussion 
and conclusions, and suggests a possible application of the present kind of 
approach to other models, most specifcally the ip\ model which is pertur- 
batively renormalizable, but generally regarded as becoming a (generalized) 
free theory when studied as the continuum limit of a conventional lattice for- 
mulation. While our proposal seems to be analytically challenging, there is 
the strong possibility that numerical Monte Carlo methods may prove useful. 



2 Overview of the Model 

Preliminaries 

The present section is devoted to a presentation of the model, while the 
following sections discuss the motivation and analysis that has led to the 
present formulation. We focus on (/pf^, n > 5, models formulated as Euclidean 
functional integrals; other models may possibly be treated by analogous pro- 
cedures. 

We suppose Euclidean spacetime is replaced by a periodic, hypercubic 
lattice with L sites on an edge, L < cx), and a uniform lattice spacing of a, 
a > 0. Let the sites be labeled by A; = {ko, ki, k2, ■ ■ ■ , ks), where kj G Z. A'o 
denotes the future time direction under a Wick rotation, and s = n — 1. We 
denote lattice sums (and products) over all lattice points by (and 11^), 
and, importantly, lattice sums (and products) over just a spatial slice at a 
fixed value of ko by E'^ (and 11'^). The total number of sites is N — L'^, while 
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the number of lattice sites in a spatial slice is A^' = L*. 

In eventually taking the continuum limit we shall do so in two steps. 
First, we let L — > oo and a — > together so that the full spacetime volume 
V = (La)"' remains large but finite; so too for the spatial volume V = {Lay. 
Second, we take the limit that both V and V diverge. In this fashion we 
can discuss finite spatial volumes which would have been less convenient if 
we had let L — oo before taking the limit a — > 0. 

Lattice action 

Following aspects of the discussion in [3], we first introduce an important set 
of dimensionless constants by 



where Sk,i is a Kronecker delta. This notation means that an equal weight of 
l/(2s + 1) is given to the 2s + 1 points in the set composed of k and its 2s 
nearest neighbors in the spatial sense only; J^^i = for all other points in that 
spatial slice. [Specifically, we define Jk,i = l/(2s + 1) for the points I = k — 
{ko, ki, k2, . . . ,ks), I = {ko, ki ± 1, k2, . . . , k^), I = (fco, /ci, /c2 ± 1, . . . , kg),. . . , 
I = {ko, ki, k2, ■ ■ . , kg ± 1).] This definition implies that T^^Jk,! = 1- 

We next write the lattice action for the full theory, including the quartic 
nonlinear interaction as well as the proposed counter term, as 



where k* denotes the n nearest neighbors to k in the positive sense, i.e., 
k* e {{ko + l,ki, . . . ,ks) T . . 1 {ko, ki, . . . ,ks + 1) }■ The last term, which 
represents the heart of the present procedure, is the suggested counter term 
and is given (with all the following sums over the spatial slice at fixed ko) by 
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^k,le{kUk„„} 1 



(1) 



2s + l 



7(0, h, N) = lEfeEfe* {^k* - <t>kf a"-2 + \mlY.,ct>l a 





Jr,kJt,k4'k 




(3) 
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Observe that we have included the proper dependence on h for the counter 
term implying that its contribution disappears in the classical limit in which 
— > 0. It may be noticed that each of the separate parts of the counter 
term scales as the inverse square of the overall field magnitude. The reason 
we have chosen the given counter term will be discussed in the following 
sections; in Sec. 6 we even show that the counter term may be considered to 
arise from a factor ordering ambiguity of the conventional theory. 

One feature of the counter term is the fact that each term involves up to 
two nearest- neighbor, spatially separated lattice points. This feature is part 
of the regularization in the lattice formulation of the model. However, if a 
second-order phase transition is achieved in the continuum limit, then such 
a regularization should still lead to a relativistic theory in the limit. 

It is important to note that as Aq — > and the quartic interaction is 
turned off, the lattice action does not pass to that of the usual free theory 
but to that of the free theory plus the original counter term. Such a theory 
has been called a pseudofree theory [4], and we shall show that the interacting 
theory with Aq > exhibits a divergence-free perturbation series about the 
pseudofree theory. A natural argument in favor of the pseudofree theory is 
given in Sec. 6. 

In conventional quantum field theory, counter terms are chosen to deal 
with the emergence of divergences; in the approach adopted in this paper, 
the counter term is chosen to deal with the cause of divergences. Relative to 
conventional treatments, therefore, it is safe to say that using the new counter 
term changes everything relative to what one normally expects based on the 
usual free theory. In particular, do not look for 'normal ordering'; instead, 
look for 'multiplicative renormalization'. 

Generating function 

One important ingredient has been left out of the lattice action, and that is 
the factor Z representing the field strength renormalization. We introduce 
this factor most simply by adopting the following expression for the lattice 
space generating function: 

S{h) = Moj e^~"^^khk(t>ka^/n - I{ct>,a,N)/h ^^^^^ ^ (4) 
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where {hk} determines an appropriate test sequence, and the normahzation 
factor Mo ensures that S{0) = 1. By a field rescahng, i.e., (pk Z^^'^cpk, the 
factor Z can be removed from the source term and introduced into the lattice 
action; we shall have occasion to use both forms of this integral. The form 
of the generating function in terms of physical fields is given in Eq. (61). 
The continuum limit will be taken as 

E{h} = lim S{h) , (5) 

a,L 

as a — > and L — oo together such that, as discussed above, C = La 
remains constant and finite. The argument of E involves suitable limiting 
test functions hk — > h{x) where ka ^ x E M". For sufficiently large C, it 
may be unnecessary to take the final limit C — > oo. 



3 Rationale for Counter Term 

Prom the lattice action it is a simple step to write down the lattice Hamilto- 
nian operator 



i^'«"^'e;^ + vo(0) + i^^E;-^.(0)' 

'fe 

+ AoEfeVt «^ + If^'Zk'W) a'-Eo. (6) 



In addition, we introduce the ground state ^(0) for this Hamiltonian op- 
erator. The constant Eq is chosen so that satisfies the Schrodinger 
equation 

^*(0) = O, (7) 
which implies that the Hamiltonian operator can also be written as 
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Since the ground state ^'(0) does not vanish, it can be written in the generic 
form 

j,-C/(,#.,a,B)/2 

and thus (using the abbreviation X,k= dX/d(f)^ and the spatial summation 
convention) 

V{4>) = \h^a-^De^l^[D-^e-^l%kj, (10) 
^\h^a-'[\U\k-\ U,kk U,k D,k +2 D--" Dl - D^' D ,kk ] ■ 

We insist that the atypical counter term J-'k{<p) is determined by the denom- 
inator D alone by requiring that 



^'EkM<l>) = W D D-\kk = a-^[2D-^D\k D,,^ ] • (11) 
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There are multiple solutions to this equation which lead to ground state 
ftmctions that are locally square integrable near the origin in field space. 
However, the only solution consistent with a nowhere vanishing ground state 
is given (up to an overall factor) by 

Dicp)^u',mj,,cpn^^'-'y'^' . (12) 

In point of fact, D was chosen first, and the counter term was then derived 
from D by this very differential equation. Why we have chosen this specific 
form for D is discussed below. 

The ground state ^'(0) leads to the probability density 

where K accounts for normalization of this expression given the additional 
assumption that U (0, a, H) — 0. The normalization integral itself then reads 

f ^-U{(t),a,h) 

Before commenting on this integral further, we wish to discuss several 
simpler integrals. 



6 



A discussion of many-dimensional integrals 

Consider the family of Gaussian integrals given (at some fixed ko) by 

Ig{2p) ^ J i^M'e-^^'k'l^l U',dcl>k , (15) 

where p G {0,1,2,3,...}, and we assume that A is of 'normal size', e.g., 
0.1 < A < 10. Although these integrals can be evaluated exphcitly, we 
prefer to study the qualitative behavior of such integrals for large values of 
N', i.e., when iV' ^ 1. For this purpose it is highly instructive to introduce 
hyper- spherical coordinates [5, 6] defined by 

(f)k = KTJk , < K <oo , -1 <r]k <l , 

J^Wk^i, K<I>1^^'- (16) 

Here k, acts as a hyper-radius field variable and the {rjk} variables constitute 
an A^'-dimensional direction field. Note well that k = ^JY!^(\)\ is the 'radius' 
of all the field variables in a given spatial slice of the lattice at some fixed ko- 
In terms of these variables, it follows that 

Ig{2p) = 2 J i^fe-^^" K'^'-'dKd^l - U'.drj, . (17) 

Observe that the integrand depends on the radius k, but it does not depend 
on the angular variables {rjk}- For very large N', the integral over k can 
be studied by steepest descent methods. To leading order, the integrand 
is narrowly peaked at values of k ~ (N' /2Ay^'^ , namely at large values of 
K. As a consequence, for each value of A, the integrand is supported on a 
disjoint set of k as A*" — >• cxd. This well-known fact [7] leads to divergences 
in perturbation calculations. For example, let us study 

^g(2) = / [K<I>1] e-^*^''^^l Kd<f>k , (18) 

which is the same type of Gaussian integral for a different value of A. For 
this study, we introduce the perturbation series 

7S(2) = 7G(2)-AA7G(4) + i(A^)2 7G(6)--- - , (19) 
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where = A* — A. Since /G(2p)//fi(2) = 0(A^'(P-i)), this series exhibits 
divergences as iV' — >• cxo. It is not difficult to convince oneself that such 
divergences are quite analogous to those that appear in quantum field theory; 
see [5] for further examples of this sort. 

If the factor k^^'"^) is removed from Iq{2p), we are led to consider 

I'ai2p) = 2 J i^'Ye-^^" d^,d{l - E',^) Il'.dv, • (20) 

Now, the integrand is broadly supported and no longer favors large k values. 
Consequently, a perturbation series evaluation of 

(2) = 2y" «:^e-^*'^' dK6{l - J:',4)Uidm 

= 7^(2) - AA7^(4) + liAAf 7^(6) - • • • (21) 

exhibits no divergences since I'q{2p) / Iq{2) = 0{N'^). 

Unlike the original integrals over integrals over the {r]k\ direction field 
variables cannot diverge under normal circumstances since each variable sat- 
isfies — 1 < % < 1. We will encounter such variables in the denominator of 
(14); however, the form of that denominator has been specifically chosen to 
ensure that such integrals converge near zero for all N' < oo. 

Relevance for the ground-state distribution 

The normalization integral for the ground-state distribution, Eq. (14), ex- 
pressed in terms of hyper-spherical coordinates, becomes 

r ^-U{Kr],a,h) 
2^ J U',mj,ir]fr'-mN' - ^'^"^D = 1 . (22) 

Note well that the absence of the factor k^^'~'^'> in this expression is a direct 
result of the counter term in the lattice Hamiltonian. which in turn gave 
rise to the denominator factor D in the ground state ^^(0). Just like the 
elementary examples in which k^^'~^^ was artificially removed, there is no 
peaking of the integrand in k due to that factor. For integrals such as 
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it is clear that the /t-dependence of the integrand is most hkely spread rather 
broadly; this conclusion would be false if the factor k,^^ had not been 
canceled by part of the term D"^. 

4 Correlation Functions and their Bounds 

In this section, following [3], we wish to show that the full spacetime corre- 
lation functions can be controlled by their sharp-time behavior along with a 
suitable choice of test sequences. 
Let the notation 

0„ = EfeMfe^jta" (24) 

denote the full spacetime summation over all lattice sites where {uk} denotes 
a suitable test sequence. We also separate out the temporal part of this sum 
in the manner 

0„ = Sfcoa0„/ = EfcoaS'^Mfe0fca* . (25) 

Observe that the notation 0^/ (with the prime) implies a summation only 
over the spatial lattice points for a fixed (and implicit) value of the temporal 
lattice value kg. 

Let the notation ((•)) denote full spacetime averages with respect to the 
field distribution determined by the lattice action, and then let us consider 
full spacetime correlation functions such as 

• • • 0«(29) ) = S (1) (2) , (29) O^'^ 0„'{2) " " " 0„'(2<j) ) , (26) 

"'0 '0 '•••'"'0 

where q > 1 and the expectation on the right-hand side is over products 
of fixed-time summed fields, for possibly different times, which are then 
summed over their separate times. All odd correlation functions are assumed 
to vanish, and furthermore, ( 1 ) = 1 in this normalized spacetime lattice field 
distribution. It is also clear that 

I 0u(2) • • • 0w(2<j) ) I < S (1) (2) (29) a^'' \ {(j)u'm 0u'(2) " " " <Pu'(^i) ) I • (27) 
At this point we turn our attention toward the spatial sums alone. 
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We appeal to straightforward inequalities of the general form 



(28) 



In particular, it follows that 



(0„,(1)0„,(2)0„,(3)0„,(4))2 < (0^,(1) 0^,(2)) (0^,(3) 0^,(4)) , 



(29) 



and, in turn, that 




(30) 



By a similar argument, it follows that 



I(0W'(1) 0«'(2) 



•••0„,(2.))|<^5l,[(05,,)]V2^ 



(31) 



which has bounded any particular mixture of spatial correlation functions at 
possibly different times, by a suitable product of higher-power expectations 
each of which involves field values ranging over a spatial level, all at a single 
fixed lattice time. By time translation invariance of the various single time 
correlation functions we can assert that 



which is defined at time kg , is actually independent of the time and, there- 
fore, the result could be calculated at any fixed time. In particular, we can 
express such correlation functions as 



Thus we see that a bound on full spacetime correlation functions may 
be given in terms of sharp-time correlation functions in the ground-state 
distribution. 



Before focusing on the limit a ^ and L — oo, let us note some important 
facts about ground-state averages of the direction field variables {%}. First, 




(32) 




(33) 



5 The Continuum Limit 
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we assume that such averages have two important symmetries: (1) averages 
of an odd number of rjk variables vanish, i.e., 

(?7fci---%2p+i) = 0, (34) 

and (2) such averages are invariant under any spacetime translation, i.e., 

{Vki ■ ■ ■ Vk2p) = iVki+i ■ ■ ■ Vk2p+i) (35) 

for any I e due to a similar translational invariance of the lattice Hamilto- 
nian. Second, we note that for any ground-state distribution, it is necessary 
that (tjI) = 1 /A^' for the simple reason that S'^ = 1. Hence, |(?7fc^/)| ^ 1/N' 
as follows from the Schwarz inequality. Since {[^'kVl]'^) = 1) it follows that 
ivlvf) = 0{1/N'^). Indeed, similar arguments show that for any ground- 
state distribution 

{Vk,---VkJ = 0{l/N'P), (36) 

which will be useful in the sequel. 

Next, we choose to study the pseudofree model, namely, when the cou- 
pling constant Aq = 0. This is as close as we can get to the free model 
itself. Unfortunately, we cannot solve Eq. (10) when V has the desired form. 
The best we can do is choose a form for U{(f),a,h) in (9) that leads to an 
approximate form of the pseudofree model. In particular, we choose 

U{<f),a,h) = (l/^)S',^,0fcAfc_;0,a2^ . (37) 

This expression is taken to be the form for U for the free model as if there 
was no counter term and consequently D was replaced by 1. This ensures us 
that the potential V(0) that follows from (10) agrees with the desired free 
model to leading order in h. Specifically, with the given choice for U and D, 
it follows that 

V(0) = |S^„0fcAfc_;A;_„0„a3^ - ^hAoN'a' + ^h^J:',J^ki(t)) a' 

+h[{N' - l)/4Ar']E',_,,^ J,,,0,/1,_„0^/[E;J,,,0,^] a^' . (38) 

We choose the matrix Ak-i so that the first term in (38) yields the desired 
gradient and mass terms in the Hamiltonian expression (6). In particular, 
we note that to match the quadratic, spatial lattice-derivative terms in the 
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Hamiltonian (typically the most singular of the quadratic terms), we can do 
so by choosing the elements of the matrix Ak^i = 0(a^*^*'+^-'). Although this 
leads to only an approximate form for the pseudofree model, it is sufficient 
for our limited purpose at present, namely, to determine the field strength 
renormalization constant Z. 



Field strength renormalization 

We now take up the question of the sharp-time averages given by 

j Z-P [E',hk<Pk afP^i^pf n',#, , (39) 

where Z denotes the field strength renormalization factor and {hk} represents 
a suitable spatial test sequence. These are exactly the kinds of expression 
that should become well behaved in the continuum limit for a proper choice 
of Z. Thus, we are led to consider 

K j z-^ [e;/^,0, af ^ nus;j,,0f](^'-^/^^' "'^^'^^ 

^2K j Z-^K'^[^',hkmaT Yi',mJkivfr'-^y^^' dKS{l-J:',r,l)U',dr]k. 

Our goal is to use this integral to determine a value for the field strength 
renormalization constant Z. To estimate this integral we first replace two 
factors with rj variables by their appropriate averages. In particular, the 
expression in the exponent is estimated by 

K^E'.^iTikAk-ima'' - K^K,iN'-'Ak-ia^' oc N' a^' a-^'+'^ , (41) 

and the expression in the integrand is estimated by 

[EihkVk - N'-^ [E',/ife af^ . (42) 

The integral over k is then estimated by first rescahng the variable — > 
/ {N' a^~^), which then leads to an overall integral-estimate proportional to 

Z-P [N'a'-^] -P N'-P [E'f^hk afP . (43) 
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Finally, for this result to be meaningful in the continuum limit, we are led to 
choose Z = N' '"^ a~^^~^\ However, Z must be dimensionless, so we introduce 
a fixed positive quantity q with dimensions of an inverse length, which allows 
us to set 



This is a fundamental and important relation in our analysis. 
Mass renormalization 

With notation where ((■)) denotes a full spacetime lattice-space average based 
on the lattice action for the pseudofree theory, an expansion of the mass term 
leads to a series of terms of the form 




(44) 




(45) 



which in turn can be expressed as 



Based on the inequality 




(47) 



valid when Aj > for all j, it follows that 



x{([SU)0^(i)aT)([S;(.)0^(.)ar)-- 



(PU)0^(.)aT>}^/^. (48) 



This leads us to consider 




- dn 



n,[EJJ,,,77f](^'-i)/2iV' 



(49) 



which, in the manner used previously, can be estimated as 




(50) 
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To make sense in the continuum limit, this leads us to identify 

ml ^ N' {qa)-^ , (51) 

with being the physical mass. Moreover, it is noteworthy that 

Zml^[N'-\qa)-^'-^^] [N' {qa)-^]m^ ^ [N' {qay]-^m^ , (52) 

which for a finite spatial volume V — N'a^ leads to a finite nonzero result 
for ZttIq. 

Coupling constant renormalization 

Wc repeat the previous calculation for an expansion of the quartic interaction 
term about the pseudofree theory. This leads us to consider terms of the form 

([AoS.^y^f), (53) 
which in turn can be expressed as 

'^o ^feW,jfc(2),...,feW HKwKw^^'] [^'fe(2)</'t(2)a1 • • • [S'^(rt a'] ) (54) 
and bounded by 

x{([S;a)</.ta)«T) ([SU)0t(^)ar) ■ • ■ {[KiA<^r)V^' ■ (55) 
This leads us to consider 

xdK6{l-J:'kril)Udrik, (56) 
which, in the manner used previously, can be estimated as 

\P ]vr/-p^sp 

([AoS-,0t«T)oc ^^;,^(,_,;j,^ , (57) 

and to make sense in the continuum limit leads us to identify 

Xo = N'^iqay-^ X , (58) 

with A being the physical coupling constant. Moreover, it is noteworthy that 

Z^Ao = [N'-\qa)-^^'-^^] [N'^iqay-^] X = [N'{qay]-^X , (59) 

which for a finite spatial volume V = N'a'^ leads to a finite nonzero result 
for Z^Xq. 
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Physical version of the generating function 

Based on the previous analysis we are led to reformulate the expression for the 
lattice space generating function (4). We first make a change of integration 
variables such that (pk — > Z^/'^cp^ leading to the expression 



where any constant Jacobian factor has been absorbed into a change of the 
overall normalization factor from Mq to M . Finally, we introduce the explicit 
form for the lattice action from (2) into (60) to yield 



This expression contains a formulation of the lattice space generating func- 
tion expressed in terms of physical fields and constants. 

Commentary 

In our final expression above there are several noteworthy points to be made. 
In a finite spatial volume V = N'a^ - which due to our hypercubic as- 
sumption for spacetime implies a finite spacetime volume V = N 0"" - the 
coefficients of the physical mass m and the physical coupling constant A are 
both finite and nonzero. We have shown earlier in this section that perturba- 
tion in both the quadratic mass term and the quartic nonlinear action leads 
to a series which is term-by-term finite when perturbed about the pseudofree 
theory. If the finite spatial volume is taken large enough (e.g., Milky Way 
sized), then a divergence-free perturbation series is established. In other 
words, the introduction of the unusual counter term has resolved any issues 
with typical ultraviolet divergences (and it is noteworthy that it has not been 
necessary to maintain an ultraviolet cutoff to achieve this goal). 

[Remark: Any theory exhibits infinite volume divergences for questions 
of a stationary nature. For example, such divergences even arise already when 




(60) 




-i [iv'^(ga)(-^)]-iE.E.. {^k^ - ^kfa-vn 



(61) 
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n = 1 and we deal with time alone, as for example with a conventional, sta- 
tionary Ornstein-Uhlenbeck (0-U) process U{t) = 2~^/'^ W{e^^), —oo < 
t < oo, where W{t), < t < oo, denotes a standard, Gaussian, Wiener pro- 
cess for which W{0) = 0, E{W{t)) = 0, and E{W {ti)W {T2)) = min(Ti,T2), 
where E denotes ensemble average. Although the 0-U paths are concentrated 
on bounded, continuous paths, it nevertheless follows that 

E(/!°^C/(t)^dt)=oo (62) 

due to the stationarity of the process.] 

Unlike the quadratic mass and quartic interaction terms, the coefficients 
of the derivative terms and the inverse-square field counter term are inverse 
to one another and do not have finite nonzero limits when L 00 and 
a — * such that the spatial volume V = [LaY is finite. This aspect is not 
unexpected since (i) the coefficient of the counter term is intimately linked 
to that of the derivative term so that the field power that appears in the 
denominator factor in the ground state distribution is {N' — 1)/2N', and 
(ii) this fact leads to a significant redistribution of probability toward the 
origin of field space, which then requires an asymptotically small Z factor to 
reestablish reasonable field averages. 



Numerical studies 

In a certain sense, the most basic pseudofrce model lacks both the quartic 
coupling and the mass term. This leads to the idealized (I) model described 

by 

Si{h) = M y"exp{E/ife0fca"/n 

-i [iV'2(ga)(-i)]-iE.E.. - Ara-yn 
-lff[N'^qaY'-'^ZkM<P)^''/^} W(f>k . (63) 

If we combine the factors of a, as well as transform the fields (pk to remove 
both h and q from the idealized action, it follows that 

Si{h)^M y"exp{E/ifc0fca" 

-iN'-'EkEk' i<l>k* - <t^kf - \N'^Y.k^km ^kd<t^k , (64) 
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where hk = g^"* ^^^"^ hk/h}^'^ and 

8k{(t>) = a''M<P) ; (65) 
note that is dimensionless apart from any that may arise from its 

inverse- square field dependence. One last transformation, in which (pk 

N' (pk, leads to the expression 

Si{h) J exp{E/ifc0fe AT'a" 

-iE.E.* - <pkr - ij:Am^kd<p, ■, (66) 

in the above expressions we have used the same symbol (M) for the nor- 
malization factor even though it has absorbed different Jacobian factors. 
This last version may be useful for numerical studies of this basic pseudofree 
model. As argued in Sec. 5, the mass term and the quartic coupling terms 
can both be added by divergence-free perturbation series. 

More directly, Monte Carlo studies may be made of the full nonlinear 
lattice theory. This can be done with a variety of rescaled field variables, but 
perhaps the most convenient is that of Eq. (66). In those field variables, the 
Euclidean lattice space probability distribution is given by 

Cexp{-iE.Efc. (0fc* - 0fc)' - k[N'aq-'WE,4>l 

-[N"aq^-']XZ,<l>t-lEMm, (67) 
where C denotes an overall normalization factor. Traditional Monte Carlo 
methods may be used with this weighting to approximately determine various 
correlation functions. The coefficients of the and A terms are direct tran- 
scriptions of those determined earlier in the present section. They have been 
estimated to be suitable for all values of a and L for which La < oo. Thus 
they should also hold in the continuum limit with finite spatial and spacetime 
volumes. Any divergences that arise in the infinite spatial and spacetime vol- 
ume limit should only be those that typically arise for stationary questions, 
such as illustrated earlier with the one-dimensional Ornstein-Uhlenbeck pro- 
cess; such divergences are expected and do not require any special treatment. 

6 Additional Discussion and Conclusions 

In [4] there is an extensive discussion of soluble, scalar nonrenormalizable 
models that are idealized versions of the relativistic model treated in the 



17 



present paper. Such models differ from the relativistic model in that in one 
case all spacetime derivative terms are omitted from the classical action and 
in the second case all but one of the derivative terms are dropped; in this 
latter case, the remaining term is identified with the eventual time direc- 
tion in an analog of a Wick rotation. These models have no physics and 
are only of academic interest. Nevertheless, from a mathematical viewpoint 
both of these models lead to Gaussian results if no further counter terms are 
introduced, and if they are studied perturbatively, they are both nonrenor- 
malizable. Fortunately, both of these models have sufficient symmetry so 
that they can be rigorously solved on the basis of self-consistency without 
using any form of perturbation theory. One of the results for both models 
is that as the coupling constant of the nonlinear interaction term is reduced 
to zero, the theories do not return to the appropriate free theory but instead 
they pass continuously to an appropriate pseudofree theory. Moreover, both 
theories exhibit meaningful, divergence-free perturbation theories about the 
pseudofree theory, but definitely not about the customary free theory. These 
models are explicitly worked out in Chaps. 9 and 10 of ref . [4] , but there is 
also a natural reason why such results are plausible. 

Nonrenormalizable quantum field theories have exceptionally strong in- 
teraction terms. This statement can be quantified as follows: Consider the 
(f^ relativistic scalar theory which has a free (Euclidean) action given by 

W ^lJ[{V(p){xY + m^ip{xY]d"x , (68) 

where x e M". These theories have the nonlinear interaction term 

V^Jip{xyd\, (69) 

where we focus on cases where p e {4, 6, 8, . . .}. Such expressions appear in 
a formal functional integral such as 

If lim;^_»o S\{h) — So{h), then the interacting theory is continuously con- 
nected to the free theory; if, instead, lim;^_^o Sx{h) — So{h) ^ So{h), then 
the interacting theory is not continuously connected to the free theory, but 
rather it is continuously connected to a psuedofree theory. Under what con- 
ditions could this latter situation arise? 
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Consider the classical Sobolev-type inequality ([4], Chap. 8) given (for 
99 ^ 0) by 

{Jip{xyrxy/P/{J[{Vip){xf + m^ip{xy]d\} < R (71) 

as a function of the parameters p and n. For p < 2n/{n — 2), it follows that 
R = 4/3; for p > 2n/{n — 2),R = 00 holds. This dichotomy is exactly that 
between perturbatively renormalizablc and perturbatively nonrenormalizable 
models. But why should this inequality relate to renormalizability? 

It is the author's long-held belief that the explanation arises from a hard- 
core behavior of nonrenormalizable interactions [8] . Simply stated, the inter- 
action for such theories is so strong that a set of nonzero measure of the field 
histories allowed by the free theory alone is projected out when the interaction 
term is present. For example, R — 00 means that there are fields for which V 
is not dominated by W in the same way as when R = 4/3. This fact suggests 
that for positive coupling constant values, some of the field histories are pro- 
jected out never to return as the coupling constant passes to zero. This result 
is the effect of a hard core at work. For the first idealized model mentioned 
above, the analogous ratio is {/ (p{x)'^ d^x}^^'^/ {/ m'^(p{x)'^ d^}, which for any 
n >1 clearly has no finite upper bound, while for the second idealized model, 
the appropriate ratio reads {J (p{xy d^x}^^"^/ {J[ip{xy+m'^ ip{x)'^] d'^x}, which, 
in this case for any n >2, has no finite upper bound. These soluble models - 
each more singular in principle than the relativistic models - are examples of 
hard-core interactions that nevertheless have divergence-free perturbations 
about their own pseuofree model. This set of facts strongly suggests that 
relativistic scalar fields such as n > 5, as we have focused on, also have 
a corresponding hard-core behavior. 

The soluble, idealized models have formulations that involve inverse field 
powers; indeed, the second model, which lies closer to the relativistic models, 
has an inverse-square field power counter term in the lattice action itself. 
For both idealized models, the needed counter term was not assumed, it was 
derived, thanks to a large symmetry of the model. For the relativistic models, 
there is insufficient symmetry to derive the needed counter term, and thus 
the counter term must be postulated, i.e, guessed. There have been several 
past suggestions that have not lived up to expectations. The present paper 
offers one more proposal that seems to satisfy the expected requirements. 

Unfortunately, the present model is not (or at least seems not) analytically 
tractable, probably lacking a technical means to analytically perform pertur- 
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bation calculations about the pseudofree theory. Nevertheless, it would seem 
possible that numerical Monte Carlo calculations should be feasible. The 
first such calculation that should be made is a test for non-triviality that 
is applied to such theories by testing whether or not the Gaussian property 
that 

((EfcMfca")') - 3((EfcMfc = (72) 

holds true for all choices of {hk} as one approaches the continuum limit; 
a single violation of this inequality would demonstrate that the continuum 
limit is not that of a free theory. In view of the connection of such full 
spacctimc correlation functions to those on a single spatial surface, as shown 
in Sec. 4, it seems unlikely that (72) holds true thanks to the chosen form of 
the counter term; on the other hand. Sec. 4 ultimately involves inequalities 
which might allow (72) to sneak through. 

Should the above nontriviality test prove successful for a (p^ relativistic 
model, for example, it would be worthwhile to study various three-dimensional 
models such as (/^g, (^g'', etc. [Remark: It is noteworthy that if instead of the 
interaction Ao0^ discussed in the present paper we had started with go4>n , 
where r G {3, 4, 5, 6, . . .}, then it follows that such an interaction would also 
possess a divergence-free perturbation expansion provided we choose 

go = N'^''-'^a'^^-'^^-'g, (73) 

where g is the physical coupling constant. Moreover, it follows that 

Z^go = [N'{qar]-'9, (74) 

for all values of r, just as was the case for the mass term in (52) and for the 
quartic coupling in (59).] It would also be interesting to reconsider the (pf 
model, which, although perturbatively renormalizable, has the property of 
tending to a Gaussian theory in the continuum limit as shown by renormal- 
ization group studies as well as Monte Carlo studies. It is possible that a 
nonperturbative and nontrivial ipf model is still to be discovered. 

The counter term from a factor ordering ambiguity 

It is not without interest that the chosen counter term can be viewed as aris- 
ing from a factor ordering ambiguity. Let tt^ denote the classical momentum 
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conjugate to the field 0^ for all A; in a spatial slice. The classical Hamiltonian 
reads 



H = Sk^l^' + ^o(0) = ^EkD^kD-^^kDa' + l^o(0) , (75) 

where D = D{(l)) = H^fEJ Jfe^^^f Passing to the quantum theory 
leads to the Hamiltonian operator [c.f., (11)] 

2 " d(t)k d(t)k 
= -f «-^EL^ + Vo(0) + ^T:,n{ct>) . (76) 

Classical limit 

When dealing with a nonrenormalizable theory, we argued in Sec. 1 
against choosing either no counter terms or those counter terms suggested 
by a regularized perturbation analysis about the free theory. This was due, 
in part, to the fact that such theories tend not to have the correct classical 
limit, namely, the original nonlinear classical theory one started with. That 
property is clear in the case of no counter terms, which leads to a free quan- 
tum theory, and it is effective as well when perturbative counter terms are 
considered because in the latter case there is no complete and well-defined 
quantum theory for which the classical limit can be studied. 

One strong test of whether or not the ideas in this paper have some 
validity would be to try to take the classical limit and confirm that the 
expected nonlinear relativistic theory emerges. The study of this question 
first requires having some control on the continuum limit, but in support of 
its possible realization we note that the second idealized model treated in 
Chap. 10 of [4] , namely, the model including the time derivative of the field, 
has been shown to indeed have the correct classical limit for the idealized 
but nonlinearly interacting model in question. It is possible that similar 
techniques may be used to show that the presently proposed quantization 
scheme for nonrenormalizable </?^ models has a classical limit that agrees 
with the original nonlinear classical theory. 

Even if our proposal leads to a nontrivial quantum theory, and even if 
that quantum theory exhibits a correct classical limit, the question may 
arise whether this proposal for quantization is the "correct" quantization 
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procedure. As in any quantization procedure, where one starts from a theory 
with ^ = and constructs a theory with h> 0, there is a great deal of latitude 
in the result. Nevertheless, in the absence of any other satisfactory proposal 
to deal with nonrenormalizable theories, one might look favorably on a model 
that offers more than was previously available. 
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